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STEADY  STATES  AND  THEIR  STABILITY  OF  HOMOGENEOUS, 
RIGID,  EXTENDED  NEMATIC  POLYMERS  UNDER  IMPOSED 
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GUANGHUA  Jit,  QI  WANG*,  PINGWEN  ZHANG§, 

HONGYUN  WANG^,  AND  HONG  ZHOUH 

Abstract.  We  study  the  steady  state  phase  behavior  of  homogeneous,  rigid,  extended  (polar) 
nematic  polymers  or  nematic  dispersions  under  imposed  magnetic  (or  electric  fields),  in  which  the 
intermolecular  dipole-dipole  and  excluded  volume  interaction  as  well  as  external  field  contribution  are 
accounted  for.  We  completely  characterize  the  phase  diagram  for  polar  nematics  with  and  without 
permanent  magnetic  moments  (or  dipoles).  For  nematics  without  magnetic  moments,  the  steady 
state  is  either  purely  nematic  or  polar  depending  on  the  strength  of  the  excluded  volume  and  the 
dipole-dipole  interaction,  in  which  the  nonzero  polarity  vector  (the  first  moment  vector)  is  coaxial 
with  the  second  moment  tensor;  thereby  the  steady  state  pdf  is  determined  essentially  by  up  to  three 
scalar  order  parameters  and  a  rotational  group  of  SO(2)  transverse  to  the  imposed  field  direction. 
For  nematics  with  permanent  magnetic  moments  (or  dipoles),  the  steady  states  are  polar  and  the 
polarity  vector  is  parallel  to  the  external  field  direction  when  a  necessary  condition  for  parameters 
is  met.  When  the  condition  is  violated,  only  stable  steady  states  have  their  polarity  vector  parallel 
to  the  external  field  direction  and  there  are  thermodynamically  unstable  nonparallel  states.  The 
stability  of  the  steady  states  is  inferred  from  the  minimum  of  the  free  energy  density. 
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bifurcation,  steady  state,  phase  transition,  polarity,  stability 
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1.  Introduction 

Kinetic  theory  is  a  useful  tool  in  modeling  soft  matter  and  complex  fluids  [2, 
4,  8,  24,  9].  In  the  past,  it  has  been  perceived  as  colossal  and  complicated — hardly 
accessible  to  theoretical  analysis.  Given  the  rising  interest  in  kinetic  theory  in  the 
mathematics  community  these  days,  various  attempts  have  been  made  to  analyze 
the  properties  of  the  partial  differential  equations  in  the  kinetic  theories  and  obtain 
their  solutions  semianalytically  and  numerically  [23,  11,  12,  13,  14,  15,  16,  17,  18,  5, 
6,  7,  10,  23,  34,  35,  26,  40,  38,  20].  A  recent  review  of  the  state  of  the  art  in  the 
mathematical  and  numerical  analysis  of  multi-scale  models  of  complex  fluids  is  given 
by  Li  and  Zhang  [25]. 
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STEADY  STATES  OF  EXTENDED  NEMATICS  POLYMERS 


In  a  series  of  papers  [38,  20,  41,  42],  we  have  given  a  systematic  study  of  the  exact 
solution  of  the  Smoluchowski  equation  with  the  Maier-Saupe  excluded  volume  poten¬ 
tial  [29,  30,  31]  and  (or)  the  dipole-dipole  interaction  potential  in  implicit  form  using 
a  reduced  order  or  projection  method  for  the  Doi-Hess  kinetic  theory  [9].  In  [38],  we 
studied  the  phase  behavior  of  rigid  nematics  dominated  by  either  the  excluded  volume 
or  the  dipole-dipole  interaction  coupled  with  an  imposed  external  field.  Therein,  we 
identified  the  order  reduction  strategy  that  we  have  adopted  in  all  the  subsequent 
studies  in  the  series.  In  [41,  42],  we  proved  a  series  of  theoretical  results  establishing 
the  conditions  for  bifurcations  of  steady  states  and  the  asymptotic  behavior  of  steady 
states  near  the  bifurcation  point  in  phase  space  in  models  with  both  dipole-dipole 
and  excluded  volume  interaction.  In  [20],  we  proved  additional  theoretical  results  for 
steady  states  of  the  Smoluchowski  equation  and  obtained  the  complete  bifurcation  di¬ 
agram  for  polar  nematics  absent  of  external  fields.  In  the  above  studies,  we  applied  the 
reduced  order  method  to  various  limiting  cases  of  dilute  or  concentrated,  nematics  or 
extended  (polar)  nematics,  with  or  without  the  external  field.  By  extended  nematics 
we  mean  here  the  nematics  that  exhibit  instantaneous  dipoles  or  magnetic  moments; 
therefore,  the  dipole-dipole  interaction  must  be  accounted  for.  We  have  established 
rigorous  mathematical  theories  for  the  existence  of  certain  types  of  solutions  and  their 
stability  [38,  20,  41,  42],  as  well  as  detailed  phase  bifurcation  diagrams  [38,  20].  One 
of  the  examples  of  extended  nematics  or  polar  nematics  is  the  nematic  liquid  crys¬ 
tal  polymer  consisting  of  bowlike  molecules,  which  were  studied  experimentally  and 
numerically  in  the  past  [3,  21,  27,  28,  39]. 

In  this  paper,  we  extend  the  previous  studies  on  homogeneous  rigid  extended 
or  polar  nematics  by  examining  the  steady  states  subject  to  the  dipole-dipole  and 
excluded  volume  interaction  under  imposed  external  fields  in  dilute  and  concentrated 
regimes  with  or  without  the  permanent  magnetic  moments  (or  dipoles).  This  will 
give  a  closure  to  the  study  of  the  steady  state  solution  behavior  governed  by  the 
Smoluchowski  equation  with  the  potential  given  by  the  following  form 


V  =  —kBT(fj.E  +  a(in})  -m  — 


3NkBT.  . 
- 2 - (mm) 


:  mm  — 


cxokBT 

2 


EE :  mm. 


(1.1) 


where  p.  is  the  strength  of  the  permanent  dipole  (or  magnetic  moment)  of  the  molecule, 
a  is  the  strength  of  the  dipole-dipole  interaction  potential  due  to  the  intrinsic  dipole 
(or  magnetic  moment)  on  the  molecule,  ao  is  the  difference  of  the  polarizability  (or 
susceptibility)  parallel  and  perpendicular  to  the  molecular  direction  known  as  the 
anisotropy,  E  is  the  external  magnetic  (or  electric)  field,  iV  is  a  dimensionless  pa¬ 
rameter  describing  the  strength  of  the  Maier-Saupe  excluded  volume  potential  and 
quadruple-quadruple  interaction,  is  the  Boltzmann  constant  and  T  is  the  absolute 
temperature  [9,  1].  In  this  model,  we  describe  the  orientation  of  the  extended  nematic 
molecule  by  a  unit  vector  m,  which  is  in  the  permanent  dipole  (or  magnetic  moment) 
direction  if  so  that  /r>0,  or  in  the  direction  of  the  averaged  dipole  when  p.  =  0 
so  that  a  >  0  in  the  latter  case.  We  remark  that  all  the  material  parameters  are  nor¬ 
malized  with  respect  to  in  this  paper.  So,  we  effectively  set  =  1  in  all  our 
numerical  calculations. 

We  adopt  the  Maier-Saupe  potential  as  the  excluded  volume  potential  in  the  stud¬ 
ies  [29,  30,  31]  because  it  is  an  approximation  to  the  integral  form  Onsager  potential 
[9].  Although  the  Maier-Saupe  and  Onsager  potential  give  the  same  behavior  qual¬ 
itatively,  theoretical  analysis  with  an  Onsager  potential  would  require  a  completely 
new  set  of  tools  which  are  yet  to  be  discovered,  because  the  Onsager  potential  re¬ 
quires  a  complete  set  of  spherical  harmonic  modes,  while  the  Maier-Saupe  only  needs 
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the  first  few.  It  is  this  property  of  the  Maier-Saupe  potential  that  results  in  our 
reduced  order  method  which  has  led  to  many  in-depth  theoretical  results  recently 
[5,  6,  7,  10,  34,  35,  38,  20,  41,  42], 

The  steady  state  solution  of  the  Smoluchowski  equation 

^^f=n-{Drfnu)  (1.2) 

in  the  Doi-Hess  kinetic  theory  with  the  normalized  chemical  potential  U  =  \nf  + 
is  given  by  [19,  9,  33,  37,  35,  38] 


=  (1.3) 

where  /  is  the  orientational  probability  density  function  (pdf)  of  the  ensemble  of 
the  extended  nematic  molecules,  Z  denotes  the  normalizing  constant  or  the  partition 
function,  V  is  the  total  potential,  77  =  mx  {d/dm)  is  the  rotational  gradient  operator, 
and  Dr  is  the  rotational  diffusivity,  which  is  taken  to  be  a  material  constant  in  this 
study. 

A  direct  consequence  from  the  above  steady  state  solution  is  that  the  total  torque, 
the  total  ensemble  averaged  torque,  as  well  as  the  averaged  external  torque  all  vanish 
at  steady  states;  i.e., 


{TZU)  =  {TZV)  =  {nVe)  =  0,  (1.4) 

where  I4  =  fcBr(— ^EE:mm  — /iE-m)  is  the  external  potential  due  to  the  external 
field  [41,  42].  In  fact,  in  steady  states, 


U  =  const. 


(1.5) 


Then, 


o=nu  =  {nu)  =  {{mnf+Tzv))  =  (tzv),  (i.e) 


since 


(771n/)  =  f  77/dm  =  0. 

7|lm|I  =  l 

In  addition,  we  denote  the  intermolecular  potential  by 

37V 

Vi  =  ksT {—a{m)  •  m - 


(1.7) 


(1.8) 


where  M=(mm).  Then, 

(77 hi)  =kBT{m  x  (— Q;(m)  —  37VM-m)) 

=  fc^T]— Q;(m)  X  (m)  —  SN (m  x  M  •  m)] 

=  —3NkBTMiiMijeijn  =  0,  (1.9) 

where  is  the  permutation  symbol  [2].  Thus, 


{nVe)  =  {'Ji{v-Vi))  =  o. 


(1.10) 
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The  steady  state  pdf  given  in  (1.3)  is  parametrized  by  the  material  and  external 
field  parameters  along  with  the  first  and  second  moment  of  m  with  respect  to  the  pdf 
defined  below, 


(m)  =  /  m/(m)(im,  (mm)  =  /  mm/(m)c?m,  (l-H) 

./|]ml!  =  l  ./l|m||=l 

which  consists  of  a  total  of  8  nonlinear,  implicit  algebraic  and  integral  equations.  In 
[20,  41,  42],  we  studied  the  case  where  the  external  field  is  absent  by  developing  an 
order-reduction  method.  We  showed  that  the  polarity  vector  (or  the  first  moment 
vector)  must  be  either  zero  or  coaxial  with  the  nematic  order  tensor  (or  the  second 
moment),  thereby  reducing  the  parameterizing  equation  of  the  Boltzmann  function 
to  three  equations  for  the  scalar  order  parameters  in  the  eigen-frame  of  the  nematic 
order  tensor.  Then,  a  comprehensive  bifurcation  diagram  for  the  equilibrium  phases 
was  obtained. 

In  the  following,  we  extend  the  order  reduction  strategy  to  extended  nematics 
under  imposed  external  fields.  We  will  present  the  phase  diagram  for  the  extended 
nematics  at  various  strengths  of  the  imposed  external  field  and  characterize  the  various 
phases  with  respect  to  the  material  and  external  field  parameters.  We  present  the 
results  first  for  the  case  without  the  permanent  dipole  or  magnetic  moment  and  then 
discuss  the  case  with  the  permanent  dipole. 

2.  Extended  nematics  without  permanent  dipoles 

When  the  permanent  dipole  is  absent,  /x  =  0.  We  will  establish  the  relation  be¬ 
tween  the  polarity  vector,  the  external  field  and  the  second  order  nematic  tensor  by 
showing  that  the  external  field  must  be  parallel  to  one  of  the  principal  axes  of  the 
nematic  order  tensor.  The  theoretical  foundation  developed  in  [20,  41,  42]  then  en¬ 
ables  us  to  extend  the  order-reduction  scheme  to  study  the  steady  state  solutions  of 
the  Smoluchowski  equation. 

First,  we  show  that  the  external  field  must  be  one  of  the  eigenvectors  of  the 
nematic  order  tensor  in  steady  state. 

Theorem  2.1.  The  external  field  is  parallel  to  one  of  the  principal  axes  of  the  nematic 
order  tensor  (i.e.,  the  second  moment  of  the  pdf). 

Proof.  It  follows  from  the  discussion  on  the  ensemble  averaged  torque  alluded  to 
in  the  introduction  that 

- (TZV)  =a(m)  x  (m)  -|-q;o(E  •  mm  x  E)  -|-3fV(m  x  M-  m) 

kbT 

=  ao(E-  mm  x  E)  -f  3iV(m  x  M  •  m)  =ao(E  •  mm  x  E)  =0,  (2.1) 


where  M=  (mm).  Then, 


j  (E-mm  X  E)jj^ 

=  11E112|1M-E|12-(M:EE)2  =  0. 

(2.2) 

This  equality  holds  if  and  only  if 

M-EIJE. 

(2.3) 

Namely,  E  is  one  of  the  principal  axes  of  the  nematic  order  tensor  M.  □ 

The  next  theorem  shows  that  the  polarity  vector  is  related  to  an  eigenvector  of 
the  nematic  order  tensor  when  it  is  nonzero. 


G.  JI,  Q.  WANG,  P.  ZHANG,  H.  WANG  AND  H.  ZHOU 


921 


Theorem  2.2.  The  polarity  vector  is  either  zero  or  must  be  parallel  to  one  of  the 
principal  axes  of  the  nematic  order  tensor. 

Proof  We  select  x,  y,  and  2  axes  such  that  the  second  moment  is  diagonal.  That 
is,  <mimj  >=0  for  i^j.  We  want  to  prove  that  (m)  is  either  zero  or  parallel  to  one 
of  the  principal  axes  of  the  second  moment  (mm) . 

First  we  denote  the  principal  axes  of  (mm)  as  ei,  02  and  63  and  'E  =  Eek  for 
l<fc<3.  Let  (m)  =riei +r2e2  +  r3e3.  Suppose  at  least  two  of  ri,  r2  and  rs  are 
non-zero  (otherwise  (m)  is  already  parallel  to  one  of  the  principal  axes  of  (mm)). 
Without  loss  of  generality,  we  assume  both  ri  >  0  and  r2  >  0  (we  can  always  rename 
the  X,  y,  z  axes  to  achieve  this).  Then  the  total  potential  is 


V{mi,m2,m3) 


-a{rimi+r2m2+r3m3)  -  {cimf  +  C2ml  +  C3ml) 
-a{rimi+r2m2+r3m3)  -  V2(mi,m2,m3), 


(2.4) 


where  Cj  =  ^{mjmj)  +  ^E^6jk.  Notice  that  V2  is  an  even  function  of  mi  and  m2. 
We  will  show  that  (mim2)  >0,  which  contradicts  the  selection  of  the  principal  axes. 


(mim2)  =  —  [  mim2exp[a{rimi+r2m2+r3m3)  +  V2(mi,m2,m3)]dr 

^  J||m||  =  l 

= -^  /  mim2exp[ar3m3  +  V2(mi,m2,m3)]x 

^  •/||m|j  — 1,  mi>0,m2>0 

{exp(Q;rimi  -|-Q;r2m2)  —  exp(— arimi  -|-Q;r2m2) 

— exp(arimi  —  ar2m2) -|-exp(— arimi  — ar2m2)}(im 

4  f 

=  —  /  mim2exp[ar3m3  +  V2imi,m2,m3)]x 

^  -'|]mj|  — 1,  mi>0,m2>0 

sinh(arimi)sinh(ar2m2)(im  >  0. 


(2.5) 


□ 

In  [38,  20,  41],  we  showed  that  the  steady  state  phase  is  nematic  for  the  dilute 
(fN  =  0)  solution  of  rigid,  extended  nematics  when  a  <  1.  Next  we  will  show  that  this 
result  remains  true  in  the  concentrated  regime  as  well. 

Theorem  2.3.  when  q;<1,  the  solution  is  nonpolar  so  that  the  external  field  is 
parallel  to  one  of  the  principal  axes  of  the  nematic  order  tensor. 

Proof.  We  assume  (m)=riei  and  denote  the  total  potential  as 


V  =  kBT(—amiri  —  V2{mi,m2,m3)).  (2-6) 

We  note  that  V2 (mi, m2, m3)  is  even  with  respect  to  all  its  variables.  We  will  show 
that  when  q;<  1,  ri  must  be  zero  so  the  only  equilibrium  is  non-polar.  We  prove  it 
by  contradiction.  Suppose  ri  >  0  (otherwise  we  can  change  the  coordinate  system  to 
achieve  this). 


ri  =  {mi) 


/g  mi  exp  [ari  mi  -I- V2  (mi ,  m2 ,  m3 )]  dS' 
/g  exp [ari mi  -I- 1^2  (mi ,  m2 , m3 )] dS” 


/swithmi>o"^i^^P[^2(”^u"^2,TO3)]sinh(arimi)d5' 

/swithmi>o®^P[^2(TOi,TO2,TO3)]cosh(arimi)dS' 


(2.7) 
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where  S  =  {m|  ||m||  =1}  denotes  the  unit  sphere.  If  a  <  0,  this  expression  yields  ri  <  0, 
contradicting  ri  >  0.  For  a  >  0  and  using  the  fact  that  tanh(a;)  <  x  for  a;  >  0,  we  have 

TOiexp[V2(TOi,TO2,TO3)]sinh(ariTOi) 

<  exp[V2(TOi,m2,TO3)]cosh(arimi)tanh(Q;riTOi)  (2-8) 

<  Q;riexp[V2(TOi,TO2,TO3)]cosh(Qfrimi). 

Substituting  this  inequality  into  the  expression  for  ri  above,  we  obtain 

ri<ari,  (2.9) 


which  is  a  contradiction  when  a  <  1  and  ri  >  0.  Therefore  ri  =  0  when  a  <  1. 

This  shows  that  the  extended  nematic  solution  is  possible  only  when  q;>1.  In 
order  to  resolve  the  steady  states  for  a  >  1,  we  denote  by  n,nj_,n*  the  three  orthonor¬ 
mal  eigenvectors  of  the  second  moment  tensor  (mm)  following  the  order-reduction 
scheme  in  [38].  We  parameterize  E,  (m)  and  m  with  respect  to  the  basis  as  follows: 

m  =  cos  On  +  sin  9  cos  j_  -|-  sin  9  sin  (pn* , 


(m)  =  Si  [cos0'n-|-sin6*'cos(()'nj_  +sin0'sin())'n*],  (2-10) 

E  =  E  [cos  0£;n-|-sin0£;  cos  ^£;n_L+sin0£;sin(/)£;n*], 

where  9' ,(j)' ,9e,4>e  are  a  couple  of  Euler  angles  for  the  first  moment  vector  and  the 
external  field,  respectively,  in  the  eigen-frame,  E  is  the  magnitude  of  the  external 
field,  and  si  is  the  order  parameter  for  polarity,  known  as  the  polar  order  parameter 
[38].  We  adopt  the  biaxial  representation  for  the  second  moment  tensor  [36,  37]: 

M  =  (mm)  =  s(nn—  ^ )  -|-/3(nj_n_L  —1/3)  -I- 1/3,  (2-11) 

o 

where  s  and  (3  are  the  two  nematic  order  parameters. 

The  steady  state  probability  density  function  is  given  by 


h  =  asi  (cos  9'  cos  9  +  sin  9'  sin  9  cos{(j)  —  (j)'))  +  ^E'^  (cos  9  e  cos  9  + 

sinO E sinO cos{4'  — (j)E)y  +  ^^[(s—  ^)(cos^0  —  1/3)  -I-  ^cos2(/)sin^0],  (2.12) 

where 

Z=  [  e^dm  (2.13) 


is  the  normalizing  coefficient  or  the  partition  function.  The  order  parameters  are 
defined  by 

Si  =  (cos0cos0'-|-sin0sin0'cos((/  — (/')), /3=  (sin^0cos2(/),s  =  (P2(cos0))  -I-/3/2, 

(2.14) 
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where  si  is  called  the  polar  order  parameter  and  P2{x)  =  ~  5  is  the  second  order 

Legendre  polynomial. 

According  to  the  above  established  results,  we  can  choose  0'  =  O,  i.e.,  identifying 
the  principal  direction  of  the  order  tensor  parallel  to  the  polarity  vector  as  n.  Then, 
9e  =  0  (e||n),  or  9e  =  ^  and  (j)E  =  0  (e||n-‘-),  |  (e||n*),  where  e  is  the  unit  vector  in 
the  direction  of  E.  These  results  effectively  reduce  the  degrees  of  freedom  for  the 
polarity  vector  and  the  nematic  order  tensor  from  8  to  3  in  the  eigenframe! 

The  stability  of  the  steady  states  is  inferred  from  the  free  energy  density  of  the 
system.  The  free  energy  density  of  the  extended  nematic  polymer  system  is  given  by 
[9] 


m=[ 

J  Itm||  =  l 


keTlnf+^+Ve 


/dm. 


From  (2.10),  we  arrive  at  the  free  energy  density  at  the  steady  state: 


m=  j 


,m||=l 

=  -kBT 
=  -kBT 


— fc^rinZ  — 


/dm 


lnZ-^M:M-|||(m)| 

lnZ-y(s^-s/3  +  /32)-|s2 


+  const. 


(2.15) 


(2.16) 


To  infer  local  stability  or  meta-stability,  we  need  to  calculate  the  second  variation 
of  the  free  energy  density  with  respect  to  both  M  and  (m)  [32].  For  global  stability, 
the  minimum  value  of  the  free  energy  density  (2)  would  be  sufficient,  though.  In  the 
stability  analysis,  the  orientational  variables  in  the  free  energy  density  should  be  (m) 
and  M,  where  the  partition  function  Z  with  h  given  by 

/i  =  a(m)  •  m-l- :  mm-l- ^EE :  mm  (2-17) 

is  accounted  for.  The  Hessian  of  the  free  energy  density  is  an  8x8  matrix.  The 
eigenvalue  of  the  Hessian  reveals  the  local  stability  properties  of  the  equilibrium  with 
respect  to  not  only  the  order  parameter  perturbation  but  also  the  directors  (the 
principal  axes  of  the  nematic  order  tensor  as  well  as  the  polarity  vector).  □ 

Reduced  Symmetry. 

Due  to  the  existence  of  molecular  dipoles,  /(m,t)  7^ /(— m,t).  However,  if 
/(m,E,t)  is  a  solution  of  the  Smoluchowski  equation,  so  is  /(U-m,U  •E,t)  provided 
U  is  an  orthogonal  matrix.  In  particular,  if  the  orthogonal  transformation  satis¬ 
fies  U-E  =  E,  the  solution  is  invariant  under  the  orthogonal  transformation.  Thus,  a 
SO(2)  symmetry  remains  when  E  7^  0,  while  SO(S)  symmetry  is  valid  only  when  E  =  0. 
The  steady  state  solution  family  of  the  Smoluchowski  equation  is  thus  parametrized 
by  the  SO(2)  rotational  group  transverse  to  the  external  field  direction  for  any  given 
pdf  parameterized  by  the  order  parameters  in  an  eigen-frame  of  the  nematic  order 
tensor. 

We  next  nondimensionalize  the  permanent  dipole  strength  and  the  anisotropy 
parameter  with  respect  to  a  characteristic  external  field  strength  Eq  as  follows: 

Xo=aoEQ,fi  =  ^j,Eo. 


(2.18) 
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When  the  permanent  dipole  is  absent,  we  choose  the  characteristic  field  strength  as 
Eq  =  E.  Then,  the  external  field  is  always  a  unit  vector  in  dimensionless  form  in 
this  case.  For  simplicity,  we  will  drop  the  tilde  on  /x  in  the  following  and  assume  it 
dimensionless. 

2.1.  Dilute  limit.  We  first  examine  the  limiting  case  N  =  0,  termed  the  dilute 
limit.  When  (m)  is  parallel  to  E,  the  polar  order  parameter  is  given  by 


Si  =  (cos0). 


where 


It  can  be  simplified  to 


/ 


^  CtSi  COS0+^  cos^  0 

z 


ze 


asiz+-^z' 


dzj 


(2.19) 


(2.20) 


(2.21) 


where  Z  =  dz.  Given  the  form  of  the  pdf  solution,  the  nematic  order 

tensor  is  uniaxial  given  by 


(mm)  =  s(nn  — 1/3)  +1/3, 


(2.22) 


where  s  =  (P2(cos0)}.  We  notice  that  the  effect  of  the  imposed  field  is  to  induce  purely 
nematic  phases  (s  ^  0)  and  facilitate  the  formation  of  the  polar  nematic  phase  (si  ^  0). 
The  left  column  in  Figure  2.1  depicts  the  phase  diagram  in  parameter  space  {a,xo)^  in 
which  the  second  order  phase  transition  curve  is  defined  by  Si  =  (cos0),l  —  a((cos^0)  — 
Si)  =  0  in  the  phase  space  {a,xo)-  The  unique  phase  is  purely  nematic  on  the  left  of 
the  curve  and  polar  nematic  on  the  right  of  the  curve.  Figure  2.1  also  depicts  the 
bifurcation  diagram  at  fixed  anisotropy  xo  =  0.5,  where  a  pitchfork  bifurcation  occurs 
at  Qfc  <  3  such  that  a  pair  of  stable  uniaxial  polar  nematic  phases  emerge  for  a>  ac- 
In  steady  states,  the  uniaxial  order  parameter  s  is  positive  for  all  a,  indicating 
that  the  material  is  in  nematic  phase  so  long  as  there  is  an  external  field  E.  Figure 
2.1  demonstrates  that  the  polar  order  indeed  enhances  the  nematic  order  [20].  The 
only  stable  solution  at  larger  a  is  the  unique,  prolate,  polar  nematic  phase  since  Si 
and  —Si  give  the  same  nematic  order  parameter  s>0. 

When  (m)  is  perpendicular  to  E,  the  polar  order  parameter  is  given  by 

si  =  ^  f  f  (2.23) 

^ Jo  J-1 

where  Z  =  J^i/q  -i)cos  0) The  pdf  solution  is  given  by 


asi  COS0+ 


2 


sin^  0cos^  (p 


(2.24) 


It  is  biaxial  if  si  ^  0,  and  the  two  order  parameters  s  and  /3  are  calculated  by 


s  — 


2 


(P2(cos0)),/3=  (sin^0cos20). 


(2.25) 


The  right  column  in  Figure  2.1  depicts  the  phase  diagram  in  the  parameter  space 
(a,Xo))  where  the  second  order  phase  transition  curve  is  defined  by  si  =  (cos0),l  — 
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(m)//E 


a 


Xo=0.5 


(m)lE 


a 


Xo=0.5 


Fig.  2.1.  The  steady  states  for  diluted  extended  nematics.  The  left  column  includes  the  phase 
diagram  and  order  parameters  when  the  polarity  vector  parallels  the  external  field,  while  the  right 
column  depicts  the  same  when  the  polarity  vector  is  perpendicular  to  the  field.  si  is  the  polar  order 
and  {s,beta)  are  the  nematic  order.  When  the  polarity  vector  is  parallel  to  the  external  field,  the 
nematic  phase  is  uniaxial  (i.e.,  fi  =  0).  Notice  that  the  nematic  order  parameters  (s,/3)  are  plotted 
only  in  the  range  of  polar  nematics,  and  the  steady  state  phase  is  nematic  so  long  as  xq^O.  The 
weak  biaxiality  in  polar  nematics  exists  in  the  perpendicular  situation. 
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q;((cos^6*)  —  s^)  =0.  The  bifurcation  diagram  for  xo  =  0-5  in  (a,si)  is  shown  in  Figure 
2.1,  indicating  that  the  effect  of  the  E  field  is  to  delay  the  formation  of  the  nematic 
phase  to  beyond  a  =  3.  Figure  2.1  also  plots  the  nematic  order  parameters  s,/3  of  the 
second  moment  tensor  in  the  regime  of  polar  nematics.  The  formation  of  the  biaxial 
phase  is  completely  due  to  the  emergence  of  the  polar  order  parameter  since  the 
solution  branch  shown  in  Figure  2.1  is  a  uniaxial  prolate  phase  (s  =  0,/3y^0)  before 
the  onset  of  the  polar  nematics.  Thermodynamically,  these  solutions  are  unstable. 
We  emphasize  that  the  nematic  order  exists  as  long  as  the  external  field  is  present 
even  when  the  polar  order  is  absent  (i.e.,  when  s  =  si  =0,/37^0).  This  is  indeed  the 
field-induced  nematic  order. 

In  summary,  the  only  stable  steady  state  under  an  imposed  external  Held  is  uni¬ 
axial,  prolate  purely  nematic  at  low  a  and  prolate  polar  nematic  at  high  a  in  which 
the  director  is  parallel  to  the  external  field  direction.  The  unstable  purely  nematic 
prolate  phase  at  low  a  and  the  biaxial  polar  nematic  one  at  high  a  exist  when  the 
external  field  is  perpendicular  to  the  polarity  vector. 

2.2.  Concentrated  limit:  steady  states  with  the  polarity  vector  parallel 
to  the  external  field.  In  this  case,  0'  =  0£;=O.  We  investigate  the  steady  states 
and  their  stability  numerically  using  a  continuation  method.  First  we  note  that  the 
governing  system  of  equations  for  (si,s,/3)  admits  the  uniaxial  solution  (si  =s  =  /3  =  0) 
when  the  external  field  is  absent.  Along  the  uniaxial  branch,  the  eigenvalue  of  the 
Hessian  in  the  second  variation  of  the  free  energy  (given  in  Appendix)  changes  sign 
at 


a  = 


3 

l-b2s' 


(2.26) 


Thus,  this  together  with  the  definition  of  si,s  defines  the  critical  strength  of  the  dipole- 
dipole  interaction  potential,  beyond  which  a  polar  nematic  branch  forms  through  a 
pitchfork  bifurcation.  That  is, 

3 

Si  =  (cos0),s  =  (P2(cos0)),a=  (2.27) 


This  case  was  studied  in  [20]  for  polar  nematics,  where  a  complete  phase  diagram  was 
obtained. 

When  a  weak  external  field  is  imposed  (xo<<l:  a  small  perturbation  of  the 
situation  studied  in  [20]),  the  phase  diagram  is  surely  perturbed.  One  prominent 
change  is  that  the  steady  states  are  all  nematic  now  showing  that  the  external  field 
eliminates  all  isotropic  states  and  forces  the  steady  states  into  nematics.  Figure  2.2 
depicts  the  phase  diagram  at  xo  =  0-01  and  xo  =  0.1,  respectively.  The  pair  of  tri- 
critical  points  connecting  the  first  order  transition  curve  to  the  second  order  one 
still  exist  except  that  they  are  perturbed  slightly.  The  triple  point  is  annihilated 
since  the  existing  phases  are  either  purely  nematic  and/or  polar  nematic.  There  are 
coexistence  regions  of  highly  aligned  purely  nematic  and  less  aligned  purely  nematic 
phases  and  purely  nematic  and  polar  nematic  phases  in  parameter  spaces.  At  high 
concentration  N  and/or  high  a  only  polar  nematic  phases  survive.  When  the  strength 
of  the  anisotropy  xo  increases  further,  even  the  tricritical  points  disappear,  leaving 
an  exclusive  second  order  phase  transition  curve  between  pure  nematics  and  polar 
nematics  analogous  to  the  two  phase  transition  curves  depicted  in  Figure  2.2  2b,  2d. 

We  now  give  a  detailed  account  of  the  steady  state  phase  behavior  shown  in 
Figure  2.2  as  we  vary  the  concentration  N  at  fixed  anisotropy  xo  =  0T  and  a  few 
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(c)  (d) 

Fig.  2.2.  The  phase  diagram  in  parameter  space  {(y,N)  when  the  polarity  vector  is  parallel  to 
the  external  field  at  xo  =  0.01  and  xo  =  0.1,  respectively,  (a)  and  (b)  depict  the  limit  of  the  metastable 
curves  and  phase  transition  curves  for  xo  =  0-01)  respectively,  (c)  and  (d)  depict  the  limit  of  the 
metastability  curves  and  phase  transition  curves  for  xo  =  0.1,  respectively.  In  this  case,  all  stable 
steady  states  are  either  uniaxial  purely  nematics  (N)  or  uniaxial  polar  nematics  (PN).  There  are 
coexistent  regions  for  nematic  and  polar  nematics  labeled  N/PN.  The  phase  transition  curve  depicted 
in  (b)  and  (d)  are  the  curve  along  which  the  energies  of  the  coexisting  phases  are  equal. 


selected  values  of  a.  Figure  2.3  depicts  the  solution  at  ;^q  =  0.1  and  a  =  1.3.  A 
stable  uniaxial  nematic  phase  with  a  small  uniaxial  order  parameter  s  exists  up  to  a 
critical  concentration  Nc^\  It  then  goes  through  a  hysteresis  bifurcation  to  yield  a 

stable  prolate  uniaxial  phase  with  a  higher  nematic  order.  At  a  higher  concentration, 

(3) 

Nc  ,  the  purely  nematic  phase  gives  way  to  a  pair  of  polar  nematic  phases  through 
a  pitchfork  bifurcation  in  si.  The  stable  polar  nematic  phase  remains  prolate  and 
uniaxial.  This  scenario  persists  to  a  critical  value  of  a.  For  instance  it  remains  at 
a  =  1.5. 
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Xo=0.1  a=1.3  Xo=0'1a=1.3 


N  N 


(a)  (b) 

Xo=0.1  a=1.3 


(c) 

Fig.  2.3.  The  phase  bifurcation  diagram  for  all  the  order  parameters  in  an  imposed  external 
field  as  functions  of  dimensionless  concentration  N .  The  parameter  values  are  a  =  1.3,xo  =  0.1.  The 
stable  solutions  when  the  polarity  vector  is  parallel  to  the  external  field  are  plotted  in  solid  curves, 
while  the  ones  when  the  polarity  vector  is  perpendicular  to  the  external  field  are  plotted  in  dashed 
curves.  The  stable  solutions  are  uniaxial  prolate  nematics  (s>0,^  =  0)  that  are  weakly  aligned  at  low 
concentration  and  highly  aligned  at  concentrations  higher  than  =4.305.  The  metastable  biaxial 
states  form  a  family  of  solutions  parameterized  by  the  rotational  group  SO(2)  at  high  concentration 
whose  major  director  is  parallel  to  the  external  field  direction.  All  others  steady  states  plotted  in 
other  patterned  curves  are  unstable,  (a)  depicts  the  steady  states  in  polar  order  parameter  si.  (b) 
depicts  the  steady  states  in  the  uniaxial  nematic  order  parameter  s  for  the  parallel  state  (solid)  and 
the  biaxial  order  parameter  (3  in  the  perpendicular  state  ( dashed ).  (c)  plots  the  steady  states  in  space 
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Xo=0.1  a=2.5  Xo=0'1  a=2.5 


N  N 


(a)  (b) 

Xn=0.1  a=2.5 


(c) 

Fig.  2.4.  The  phase  bifurcation  diagram  for  all  the  order  parameters  in  an  imposed  external  field 
as  functions  of  N.  The  parameter  values  are  a  =  2.5, xo  =  O-l-  '^he  stable  solutions  when  the  polarity 
vector  is  parallel  to  the  external  field  are  plotted  in  solid  curves  while  the  other  stable  solutions  are 
plotted  in  dashed  curves.  The  curves  of  the  other  patterns  are  for  unstable  solutions.  The  stable 
solutions  are  uniaxial  prolate  ones  (s>0,/3  =  0j  which  are  weakly  aligned  at  low  concentration.  A 
stable,  uniaxial,  polar  nematic  steady  state  emerges  beyond  a  critical  concentration  Nc  as  the  result 
of  a  saddle  node  bifurcation  in  the  phase  space  (N,si),  with  the  polarity  vector  parallel  to  the  external 
field  direction.  At  a  slightly  higher  critical  concentration,  a  metastable  biaxial  steady  state  emerges 
also  as  a  result  of  saddle  node  bifurcation,  whose  major  director  is  parallel  to  the  external  field 
direction.  Other  steady  states,  uniaxial  and  biaxial  alike,  are  unstable,  (a)  depicts  the  steady  states 
in  polar  order  parameter  si.  (b)  plots  the  steady  states  in  the  uniaxial  order  parameter  s  in  parallel 
state  (solid)  and  the  biaxial  order  parameter  0  in  the  perpendicular  state  (dashed),  (c)  plots  the 
steady  states  in  the  phase  space  (N,s/ 0,si). 
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Xq=0.1  a=3.5  Xq=0'1  a=3'5 


N  N 


(a) 


(b) 


Xo=0.1  a=3.5 


(c) 

Fig.  2.5.  The  phase  bifurcation  diagram  for  all  the  order  parameters  in  an  imposed  external 
field.  The  stable  solutions  are  plotted  in  solid  (polarity  parallel  to  the  external  field)  and  dashed 
(polarity  perpendicular  to  the  external  field)  curves.  The  parameter  values  are  q  =  3.5,xo  =  0.1. 
The  globally  stable  solutions  (solid  curves)  are  uniaxial  prolate  ones  (s>0,P  =  0)  which  are  highly 
aligned.  The  metastable  steady  states  (dashed  curves)  are  biaxial  with  the  major  director  parallel  to 
the  external  field  direction.  The  others  are  all  unstable  and  are  plotted  in  curves  other  than  solid 
and  dashed  curves,  (a)  depicts  the  steady  states  in  polar  order  parameter  si.  (b)  plots  the  steady 
states  in  the  uniaxial  order  parameter  s  in  parallel  state  (solid)  and  the  biaxial  order  parameter  0 
in  the  perpendicular  state  (dashed),  (c)  plots  the  steady  states  in  the  phase  space  (A^,s//3,si). 


At  higher  strength  of  a,  say  a  =  2.5,  a  different  scenario  emerges.  The  stable 
uniaxial  nematic  phase  gives  away  to  a  pair  of  unstable  polar  nematic  phases  at  JVc^\ 
and  then  regains  stability  via  saddle  node  bifurcations  to  a  pair  of  polar  nematics  at 
lower  concentrations  shown  in  Figure  2.4.  The  stable  polar  nematic  branch  and  the 
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a=1.5xQ=0.5 


a=1.5Xp=0.5 


N  N 


(a)  (b) 

a=1.5Xg=0.5 


(c) 

Fig.  2.6.  The  phase  bifurcation  diagram  for  the  order  parameters  in  an  imposed  external  field. 
The  parameter  values  are  a  =  1.5,xo  =  0.5.  The  stable  steady  states  with  the  polarity  parallel  to  the 
external  field  are  plotted  in  solid  curves,  while  the  metastable  steady  states  with  the  polarity  vector 
perpendicular  to  the  external  field  are  plotted  in  dashed  ones.  The  globally  stable  steady  states  are 
uniaxial,  prolate.  The  metastable  ones  are  biaxial  with  the  major  director  parallel  to  the  external  field 
direction.  The  biaxial  nematic  branch  is  first  born  purely  nematic  at  an  intermediate  range  of  the 
nematic  concentration  N .  Soon  after  it  is  born,  a  pitchfork  bifurcation  in  the  polar  order  parameter 
leads  to  a  polar  nematic  phase.  The  other  steady  states,  uniaxial  and  biaxial  alike,  are  unstable,  (a) 
depicts  the  steady  states  in  the  polar  order  parameter  si.  (b)  plots  the  steady  states  in  the  uniaxial 
order  parameter  s  in  parallel  state  (solid)  and  the  biaxial  order  parameter  0  in  the  perpendicular 
state  (dashed),  (c)  plots  the  steady  states  in  the  phase  space  (N,s/0,si).  The  increase  in  xo  results 
in  the  delay  of  the  formation  of  metastable  biaxial  polar  nematics  to  the  higher  concentration. 
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a=2.8  V  =0.5  a=2.8  x  =0.5 


(c) 

Fig.  2.7.  The  phase  bifurcation  diagram  for  all  the  order  parameters  in  an  imposed  field.  The 
parameter  values  are  a  =  2.8,xo  =  0.5.  The  stable  solutions  are  plotted  in  solid  (parallel  polarity) 
and  dashed  (perpendicular  polarity)  curves.  The  stable  steady  states  in  the  case  where  the  polarity  is 
parallel  to  the  external  field  are  uniaxial  prolate  (s>0,f3  =  0)  and  highly  aligned.  At  high  concentra¬ 
tion,  a  metastable  biaxial  steady  state  emerges  with  the  polarity  vector  perpendicular  to  the  external 
field.  All  other  curves  represent  unstable  steady  states,  (a)  depicts  the  steady  states  in  polar  order 
parameter  si.  (b)  plots  the  steady  states  in  the  uniaxial  order  parameter  s  in  parallel  state  (solid) 
and  the  biaxial  order  parameter  (3  in  the  perpendicular  state  (dashed),  (c)  plots  the  steady  states  in 
the  phase  space  {N,s/ (3,si) . 


unstable  one  connected  to  it  remain  prolate  and  uniaxial.  This  persists  through  a  =  3. 
For  a  >3,  the  prolate,  uniaxial  polar  nematic  exists  beginning  dA  N  =  0  and  becomes 
the  only  stable  state.  See  Figure  2.5.  We  remark  that  the  polar  nematic  steady  states 
exist  in  pairs  with  the  polar  order  zb^i.  The  solution  behavior  changes  slightly  when 
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a=1.5)'j=1.0 


N  N 


(a) 


(b) 


a=1. 5x^=5 


1^=5  a=1 .5 


N  N 


(c)  (d) 

Fig.  2.8.  The  polar  order  parameter  si  as  well  as  the  nematic  order  parameters  sjfd  at  q;  =  1.5 
and  xo  =  l)5.  (a)  and  (b)  depict  the  order  parameters  at  xo  =  ^  (d)  depict  the  order 

parameters  at  xo  =  5-  The  effect  of  xo  steady  state  with  parallel  polarity  is  to  promote  the 

formation  of  the  polar  order  at  low  concentration,  whereas  the  effect  of  xo  state  with  perpen¬ 

dicular  polarity  is  to  delay  the  formation  of  the  metastable  nematic  order  to  higher  concentration. 


Xo  increases  to  xo  =  0-5.  The  hysteresis  bifurcation  along  the  purely  nematic  branch 
disappears.  However,  the  pitchfork,  a  transcritical  bifurcation,  persists,  yielding  the 
pair  of  polar  nematics.  Figure  2.6  and  2.7  depict  the  solution  behavior  at  a  =  1.5, 2.8. 
When  a  >3,  the  behavior  is  analogous  to  what  is  depicted  in  Figure  2.5.  When  xo 
is  increased  further,  the  general  trend  is  that  the  nematic  order  increases  and  the 
critical  concentration  and  strength  of  dipole  potential  decreases.  Figure  2.8  depicts 
two  steady  states  at  a  =  1.5  and  xo  =  1,5,  respectively.  At  the  larger  value  of  xo,  tbe 
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critical  concentration  where  the  polar  order  emerges  decreases  below  N  =  2.  The  other 
features  alluded  to  for  smaller  xo  remain. 


Io=0.1 


(a) 


8 

6 


Z4 

2 


*  Purely  Nematics  and 
l^erpendicular  Polar  Nematics 

\ 

\ 

\ 

\ 


Perpendicular 

Nematics 


Parallel  Nematics 


1  1.5  2  2.5  3  3.5 


(b) 


a 


Fig.  2.9.  The  phase  diagram  for  extended  nematics  in  an  imposed  magnetic  field  at  xo  =  0-l. 
Both  steady  states  of  polarity  parallel  and  perpendicular  to  the  external  field  are  included,  (a)  depicts 
the  limits  of  the  metastable  curves,  while  (b)  depicts  the  limit  curves  for  the  steady  states  with  the 
perpendicular  polarity.  Parallel  nematics  mean  the  uniaxial  steady  states  with  the  polarity  parallel  to 
the  external  field,  while  the  perpendicular  nematics  imply  the  biaxial  steady  states  with  the  polarity 
perpendicular  to  the  external  field.  The  curves  plotted  in  (b)  are  part  of  those  in  (a).  This  gives  all 
the  possible  phases  in  parameter  space  {a,N). 


2.3.  Concentrated  limit:  steady  states  with  the  polarity  vector  per¬ 
pendicular  to  the  external  field.  From  the  limiting  case  A^  =  0,  we  learned 
that  there  exist  biaxial  steady  states  when  the  polarity  vector  is  perpendicular  to  the 
imposed  external  field  direction;  however,  these  are  unstable.  In  [38],  we  showed  that 
stable  biaxial  nematic  steady  states  can  exist  when  excluded  volume  and  external  field 
potential  are  coupled.  In  the  dipole-dipole  interaction,  the  excluded  volume  interac¬ 
tion  and  the  external  field  are  all  coupled,  we  show  numerically  that  the  biaxial  state 
can  be  stabilized  in  the  non-dilute  limit  by  the  excluded  volume  potential,  but  these 
steady  states  are  only  metastable. 

In  the  weak  external  field.  Figure  2.3  shows  the  phase  transition  of  a  stable 
purely  nematic  biaxial  steady  state  to  a  pair  of  biaxial,  polar  nematic  phases  through 
a  pitchfork  bifurcation  in  the  polar  order  parameter  given  by  the  dashed  curves. 
The  stable  biaxial  phase  is  borne  out  of  a  saddle  node  bifurcation  with  one  branch 
metastable  at  intermediate  concentration  N.  The  existence  of  the  metastable  purely 
nematic  biaxial  state  is  reminiscent  of  the  similar  steady  state  reported  in  [38]  without 
the  dipole-dipole  interaction.  We  note,  however,  that  the  stable  biaxial  steady  state 
is  nearly  uniaxial  with  the  major  director  in  the  direction  of  the  external  field!  This 
scenario  persists  to  a  critical  value  of  a.  In  Figure  2.4  (a  =  2.5),  the  biaxial  steady  state 
is  born  stable  and  polar  through  a  saddle  node  bifurcation  at  lower  concentrations.  For 
a  >  3,  the  stable  biaxial  polar  nematic  phase  exists  for  all  fV  >  0.  The  dashed  curves 


G.  JI,  Q.  WANG,  P.  ZHANG,  H.  WANG  AND  H.  ZHOU 


935 


in  these  figures  represent  the  stable  steady  states  corresponding  to  the  perpendicular 
external  field.  So  the  effect  of  a  perpendicular  external  field  is  to  induce  metastable, 
biaxial  steady  states  with  their  major  director  aligned  in  the  external  field  direction. 
As  the  strength  of  the  external  field  increases,  i.e.,  as  xo  increases,  the  formation  of  the 
polar  nematic  phase  is  pushed  to  a  higher  concentration  at  fixed  a  and  the  metastable 
biaxial  nematic  phase  tends  to  be  born  polar  in  this  range  of  xo-  See  Figure  2.8  for 
the  case  of  a=1.5  and  xo  =  l)5.  So,  in  high  external  fields,  the  metastable  biaxial, 
steady  state  survives  only  in  high-concentration  regions  as  polar  nematic.  The  phase 
diagram  in  parameter  space  {a,N)  for  the  perpendicular  external  field  at  xo=0.1  is 
given  in  Figure  2.9,  where  the  dotted  curves  define  the  limits  separating  the  biaxial 
polar  nematics  from  the  other  nematics.  We  note  the  existence  of  a  coexistence  region 
where  pure  nematics  and  perpendicular  polar  nematics  coexist.  The  region  labeled 
parallel  nematics  is  the  one  where  nematic  solutions  exist  only  when  the  polarity 
vector  is  parallel  to  the  external  field. 

We  next  investigate  the  extended  nematics  with  permanent  dipoles  (or  magnetic 
moments). 

3.  Extended  Nematics  with  permanent  dipoles  (/xt^O) 

We  first  establish  a  lemma  to  reveal  the  relation  among  the  nematic  order  tensor, 
the  polarity  vector  and  the  external  field  in  this  case. 

Lemma  3.1.  In  steady  states,  there  exists  a  constant  c  such  that  /i(m) -|-xoM-E  = 

cE. 

Proof.  From  the  introduction, 


(We)  =  0, 

(3.1) 

so  we  have 

/x(m)  X  E-I-xoM-E  X  E  =  0. 

(3.2) 

This  implies 

/x(m)-|-xoM-E  E. 

(3.3) 

Namely,  there  exists  a  constant  c  such  that 

/x(m) -l-xoM  •  E  =  cE. 

(3.4) 

With  this  lemma,  we  can  establish  the  following  theorem. 

Theorem  3.1.  When  /xy^O,  the  polarity  vector  is  nonzero  in  steady  states.  Thereby, 
the  steady  states  are  all  polar. 

Proof.  We  show  this  by  a  counterproof.  Suppose  (m)  =  0;  the  lemma  implies  E 
is  an  eigenvector  of  the  second  moment.  We  choose  the  coordinates  such  that  the 
second  moment  is  diagonal  and  E  =  i?(0,0,l).  We  denote 

^(^)=/||m||=i"^3e'‘^^^dm, 

3 

h{X)  =  XyiEms  + 


(3.5) 
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Then  (m3)=F(l)  =  0  and  F'(A)>0.  It  can  be  easily  verified  that  T’(0)  =  0.  Thus, 
there  exists  a  0  <  A*  <  1  such  that 


F'{X*)=nE  [  mle'^dm  =  0.  (3.6) 

This  contradicts  F'{X)  >0.  The  polarity  vector  must  be  nonzero.  □ 

Analogous  to  the  results  established  in  the  previous  section,  we  can  show  the  next 
result. 


Theorem  3.2.  When  q;<1,  the  nonzero  polarity  vector  and  the  external  field  must 
be  parallel  to  each  other. 

Proof.  We  assume  E  =  i?(0,0,l)  and  (m)  =  (ri,0,r3)  with  ri>0.  Then, 


/ 


1 

Z 


cxrimi+V2(mi,m2,m3) 
^  5 


(3.7) 


where  V (mi, m2, ms)  is  even  in  mi. 


,  ,  omisinh(arimi)exp[y2(mi,m2,m3)]dm 

=  (ttii)  =  — - =  — - - - - ^ - 

/||m||=i/^"^  /|j^ll=i_„^>ocosh(arimi)exp[y2(mi,m2,m3)]dm 

^  /||m||=l,mi>0’^l‘^O®^(“’'l"^l)®^P[^2(wi,m2,m3)]dm  ^ 

/||m||  =  l.mi>0‘^O®^(“^l’^l)®^P[^2(TOl,TO2,W3)]c?m 


When  q;<  1,  ri  =0.  Namely,  the  polarity  vector  is  parallel  to  the  external  field. 

A  simple  argument  given  below  rules  out  the  possibility  of  the  existence  of  a  steady 
state  whose  polarity  vector  is  perpendicular  to  the  external  field.  Let’s  assume  r^  = 
(m3)  =  0  in  the  steady  state.  Then,  a  simple  calculation  shows 


»’3  =  (w3)  =  4  /  m3  exp 

JlimlNl 


+  ar  1  mi  +  ^  ^  M^mf  +  ^E'^m 


fxms 


2  I  Xo  r,2  2 


dm 


2 

Z 


/  m3sinh(/im3)exp 

3 

1  -  2  1  Xo  p2  2 

arimi  +  ^  >  Mum^  -I-  -—E  m3 

'  ||m||  =  l,m3>0 

z  z 

L  J 

dmfi^O. 

(3.9) 


Thus,  the  polarity  vector  will  never  be  perpendicular  to  the  external  field  in  steady 
states.  □ 

In  addition,  we  can  show  that  the  polarity  vector  and  the  external  field  are  parallel 
in  certain  parameter  regimes  given  below. 

Theorem  3.3.  When  |^|  >XoE  and  3Np.<  \axoE\,  the  polarity  vector  and  the  elec¬ 
tric  field  are  parallel  to  each  other  and  are  both  eigenvectors  of  the  nematic  order 
tensor. 

Proof  We  select  the  z-axis  to  be  the  direction  of  E  and  the  y-axis  to  be  per¬ 
pendicular  to  the  plane  spanned  by  (m)  and  E.  If  (m)  and  E  are  parallel  to  each 
other,  they  must  be  one  of  the  eigenvectors  of  the  second  moment  by  the  lemma.  We 
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therefore  focus  on  the  case  where  they  are  not  parallel  to  each  other.  In  this  Cartesian 
coordinate  system,  we  have 


m  =  (mi,m2,TO3),  I]  =  E{0,0,1),  (m)  =  (ri,0,r3), 
[/(m)  =  —an mi  —  {fiE  +  ar3)m3  —  :  mm  — 

/(m)  =  ^exp 


Z  =  I  exp 


arimi  +  {^j,E  +  ar3)m3  +  :  mm+  ^  E"^  m^ 

arimi  +  (^i?  +  ar3)m3  +  :  mm+  ^  i^^mg 


dm.  (3.10) 


Note  that  our  choice  of  the  Cartesian  coordinate  system  makes  E  nonnegative  and 
M23  =  0  from  the  lemma.  We  also  select  the  positive  direction  of  the  polarity  vector 
so  that  /i  is  always  nonnegative.  The  nonlinear  integral  equations  for  ri  and  rs  are 


ri  =  (mi)=  /  mif{m)dm, 

Js 

r3  =  {m3)=  /  m3f{m)dm. 
JS 


(3.11) 


We  use  the  spherical  coordinate  system  where  the  pole  is  the  y-axis  and  we  denote  it 
by  (V’)C)-  We  have 

(mi, m2, m3)  =  (sin'0sinC,cos'i/',sin'!/)COsC), 

U{'4>X)  =  —oiTi  sini/'sin^  —  (/xi?  +  ar3)  sinf/^cosC—  ^  sin^i/>cos^C 

SN 

— —  ( Ml  1  sin^  i/'  sin^  C  +  -^22  cos^  f/'  +  -^33  sin^  ip  cos^  ( 


+2Mi2  sinf/)  cosf/' sin  C  +  Mi3  sin^  •(/>  sin  2C) 

=  C^Mut  (■0.  C)  +  C^Ext  (■*/'>  C)  > 


(3.12) 


where  the  mutual  interaction  and  the  external  interaction  parts  of  the  potential  are 

C7Mut(V',C)  =  —Ctrl  sin0sinC  —  ar3  sini/)cos0 
3iV 

— —  ( Ml  1  sin^  sin^  ^  +  M22  cos^  0  +  M33  sin^  0  cos^  C 


+2Mi2sin0cos0sinC  +  Mi3sin^0sin2C), 
CExt(0,C)  =  — M^sin0cosC—  ^E"^  sin^0cos^C. 


(3.13) 


Ti  and  r3  are  given  by 


p7T  p‘2.'K 

ri  =<  siiiT/jsin^  >=  /  /  sin7/^sin(^/(^,^)(i?/jsin^(iC, 

JO  JO 

p-K  p'llT 

r3  =<sin0cosC>=  /  /  sin0cos0/(0,0)d0sin0c?0  (3-14) 

Jo  Jo 

We  need  to  prove  that  all  solutions  of  (3.11)  satisfy  ri  =0.  We  prove  it  by  contra¬ 
diction.  Suppose  there  is  a  solution  of  (3.11)  satisfying  n  00.  It  is  straightforward  to 
show  that  if  (ri,r3)  is  a  solution  of  (3.11),  then  (— ri,r3)  is  also  a  solution  of  (3.11). 
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Without  loss  of  generality,  we  assume  ri  >  0.  We  are  going  to  show  that  fj,  >  xoE  and 
d 

ri>0  lead  to  (^C/Ext(V’)C))  > 0;  which  contradicts  (^)  =  0  proved  in  [42].  In  fact 


(^)=0  implies  (^^)  =  0. 

Noticing  that  C/Ext(V')C)  is  an  even  function  of  C  and  consequently 
is  an  odd  function  of  C,  we  have 


C^Ext(V’>C) 


(^%xt(V',C))  =  y  J  ^C^Ext(V’,C){exp[-C/(V',C)]-exp[-C/(f/',-C)]}c^CsinV'dV'- 

(3.15) 


In  the  above,  the  first  factor  of  the  integrand  satisfies 

d 

^C^Ext('0,C)=M^^sin'i/'sinC  +  Xo£^^sin2V'cosCsinC 

dC 

=  {fi  +  xoEsinipcosQEsmtpsm^  ,  for^  G  (0,7r),  ip  G  (0,7r).  (3.16) 
The  second  factor  of  the  integrand  satisfies,  for  ^  G  (0,7r),  G  (0,7r), 

exp[-C/(V',  0]  -  exp[-C/  {-ip,  -C)] 


3 

=  2exp[^^E‘^sin^'tpcos'^(+{^J,E  +  ar3)sm'tpcos(+  ■ 

i=l 

sinh(Q;ri  sin-0  sinC  +  3iVMi2sini/)  cos  V'sinC  +  ^Mi3sin^0sin2C) 

Introducing  the  change  of  variable  ■0^7r  — 0  while  0G  [7r/2,7r],  the  integral  reduces 
to 

(^C^Ext(0,C))  =4/o’'^Vj' ^^Ext(0,C)exp[f^£;^sinVcos2C+ 

3 

{fiE  +  ars)  sin'i/^cos^  +  .  (3.17) 

sinh(Q;ri  sin0sinC  +  ^Mi3sin^'ipsin2()cosh{3NMi2Smtpcostpsm()d'ipd(. 

From  the  lemma,  we  have 

dri  =  -XoEMi3.  (3.18) 


ari  sin0sinC  H - Mi3sin^0sin2C  =  Mi3sin0sinC  (  —  *^^ill  +  3iVsin0cosC  )  . 

2  \  n  J 

(3.19) 

When  ^>|xo7?|  and  3N/r  <  |q;xo711|,  the  above  integral  cannot  be  zero  unless  ri  = 
Mi3  =  0.  Therefore,  the  direction  of  (m)  must  be  parallel  to  that  of  E.  When  these 
two  vectors  are  parallel  to  each  other,  the  lemma  implies  that  they  must  be  one  of 
the  eigenvectors  of  the  nematic  order  tensor.  □ 
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Remark  3.4.  The  inequalities  in  the  above  theorem  also  imply  that  a>iiN. 

When  the  condition  in  the  theorem  is  violated,  nonparallel  solutions  may  exist. 
Without  loss  of  generality,  we  set  E  =  i?(0,0,l),  (m)  =  (ri,0,r3),  and  parameterize 

m=  (sindcos(/),sin6*sin^,cos0).  (3.20) 

Then, 

/=ie^ 

h  =  fiE  cos  d  +  a  (r  1  sin  d  cos  (p  +  rs  cos  6)  +  ^  [Mu  sin^  9  cos^  (f>  +  Mi2  sin^  9sm2(j)+ 
Mi3  sin 29 cos (p  +  M22  sin^  9 sin^  (p  +  M33  cos^  9]  +  cos^  9,  (3-21) 

ri  =  (sin  9  cos  (p),r^  =  (cos  9) ,  Mi  1  =  (sin^  9  cos^  cp) ,  M12  =  |  (sin^  9  sin  2^) , 


Mi3  =  I  (sin  29  cos  (p) ,  M22  =  (sin^  9  sin^  (p) ,  M33  =  (cos^  9) . 

The  pdf  steady  state  solution  can  further  be  simplified  in  the  case  of  ri  7^0.  The 
following  lemma  shows  that  M12  =  0  if  ri  7^  0  indicating  that  nonparallel  steady  states 
are  in  plane  with  the  polarity  vector  and  the  external  field. 

Lemma  3.2.  When  {m2)  =  M2d,  =  Q,  (mi)Mi2=0. 

Proof.  We  note  that  the  Cartesian  coordinate  and  the  spherical  coordinate  are 
related  by  the  following  relations: 

Xi  =/3sindcos(/),  X2  =psindsin0,  X3  =  pcosd, 


^=sin6»cos^,  ^=sin6»sin^,  J(^=cos6», 

|^  =  yCOs6lcos((i,  ■^  =  ^cos9sm(p,  ■^  =  -]^sm9, 


(3.22) 


It  follows  that 


d(f)  _ 

sintp 

d4> 

_  cos^ 

d4>  _ 

dx\ 

psin9  ’ 

dx2 

ps'md  ’ 

8x3  ~ 

df 

p27r 

din=  /  dcj) 

fl 

^df  d9 

-_i  dm2 

Jo 

10 

d9  dm2 

d(p  ( 

^27r 


d(p 


df\ 

siii9cos9sm<p^ +cos(p—  ) 
oO  d(p ) 


^df 

iutiincf 

c27I- 


d9 


=  2{m2)+  f  d9  f  "^(/cJ^)^^  =  2(^2). 
Jo  Jo  d(p 


(3.23) 


If  (m2)  =  0,  {m2mf)  =  0,  the  pdf  has  the  following  form: 

ft,  =  a(rimi  +r3m3)+/xm3  +  ^mg 
3N 

+  —  [Miimf  +  M22ml  +  (I  -  Mn  -  M22)m\  +  2Mi2mim2  +  2Mi3mim3] . 


(3.24) 
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Then, 


df 

- —  =3N{M22m2  +  Mi2mi)f. 
dm2 


So  we  have 

3iVM22(m2)  +  3iVMi2(mi)  =  2(m2),  (3.25) 

i.e.,  riMi2  =  0.  □ 

We  next  examine  the  possible  nonparallel  steady  states  and  their  stability  in  the 
dilute  and  concentrated  limit,  respectively. 


Fig.  3.1.  The  critical  curves  of  a*  versus  fi 
state  polarity  vector  parallels  the  external  field, 
of  the  curves. 


Xo=1.0 


/or  xo  =  0.2,1.  On  the  left  of  the  curves,  the  steady 
Nonparallel  steady  states  may  emerge  on  the  right 


3.1.  Dilute  limit.  When  N  =  0,  the  pdf  solution  is  parameterized  by  r'i,r3, 
which  are  determined  by  the  solutions  of 

ri  =  (sin0cos(/)),r3  =  (cos0).  (3.26) 

The  parallel  versus  nonparallel  scenario  corresponds  to  ri  =  0  vs  ri  ^  0.  This  is  studied 
exhaustively  in  [42] .  It  is  shown  that  the  polarity  vector  is  parallel  to  the  external  field 
when  fi  >  Ixorl.  In  addition,  it  is  shown  that  there  exists  an  a*  such  that  when  a>  a* 
and  n  <  1x0^1  a  nonparallel  solution  exists.  Of  course,  the  condition  a{xoE  +  >  0 

must  be  met  if  the  nonparallel  solution  has  ri  >  0. 

The  bifurcation  of  ri  away  from  zero  is  accomplished  through  a  pitchfork  bifurca¬ 
tion  in  ri  as  a  >  1  varies.  The  Jacobian  of  the  two  nonlinear  equations  defining  ri,r2 
is  given  by 


J  = 


1  —  a(sin^0cos^(())  — a(sindcos0cos(()) 
— a(sindcos0cos(())  1  — a(cos^0) 


(3.27) 
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When  ri  =  0,  (sin6*cos0cos(/))  =  0.  At  the  bifurcating  point,  1  —  a(sin^0cos^(/))  =  0. 
This  defines  the  critical  strength  a*  where  the  polarity  vector  begins  to  be  non¬ 
parallel.  Figure  3.1  depicts  a  set  of  critical  curves  for  selected  values  of  xo  =  0.2  and 
1,  respectively.  The  critical  curve  in  (/r,a*)  shifts  to  higher  values  in  a*  when  |xo 
increases. 


Xj=1.0n=0-5a=2,5  N=4.75 

1i - ^ ^ ^ 


0.5 


0 


-0.5 


^ ^ ^ - 

-2-1012 

E 


X„=1.0n=0-5a=2.5N=4.75 


3Cg=1.0  p=0.5  a=2.5  N=4.75 


3Cg=1.0  p=0.5  a=2.5  N=4.75 


Fig.  3.2.  The  polarity  vector  as  well  as  the  nematic  order  tensor  components  in  concentrated 
regimes.  The  parameter  values  are  xo  =  l»c^  =  2.5,A^  =  4.75,/i  =  0.5.  The  stable  and  metastable  steady 
states  are  uniaxial,  prolate,  polar  nematics  with  the  polarity  vector  parallel  to  the  external  field.  A 
metastable  state  can  exist  at  intermediate  range  of  the  field  strength  with  the  polarity  axis  opposite 
to  the  external  field  direction.  The  nonparallel  steady  states  in  which  the  polarity  is  not  parallel  to 
the  external  field  can  exist,  but  are  unstable.  The  nematic  order  of  the  globally  stable  steady  state  is 

higher  than  that  of  the  metastable  one.  Here  ||M  —  1 1|2  =  kly^*  stable  and  metastable  states  are 
drawn  in  solid  curves.  Nonzero  values  in  ri  and  M13  indicate  the  nonparallel  nature  of  the  states. 


3.2.  Concentrated  limit.  When  A^>0,  the  existence  of  the  nonparallel 
steady  states  becomes  complex  when  the  parameters  violate  the  conditions  in  the 
theorem.  We  look  into  this  numerically.  Our  numerical  results  show  that  nonparallel 
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steady  states  may  exist;  however,  they  are  all  unstable  and  there  may  exist  up  to  two 
stable  steady  states  exhibiting  uniaxial  symmetry. 

Figure  3.2  depicts  a  representative  solution  diagram  focusing  on  the  stable  steady 
states  and  on  some  nonparallel  ones  at  TV  =  4.75.  Two  stable  uniaxial  polar  nematics 
may  coexist  in  a  finite  range  of  E.  One  exists  only  at  weak  field  strength  while 
the  other  exists  at  all  strengths.  The  short-lived  stable  state  is  metastable  and  its 
polarity  vector  is  opposite  to  the  direction  of  the  imposed  external  field,  in  contrast 
to  the  globally  stable  state.  In  all  the  numerical  experiments  we  have  conducted, 
the  nematic  order  of  the  stable  state  is  always  higher  than  the  metastable  state. 
Non-parallel  steady  states  can  exist  in  certain  parameter  ranges.  However,  they  are 
unstable!  N  =  4.75  is  within  the  biphasic  region  in  the  case  of  pure  nematics.  At  higher 
concentration,  e.g.,  N  =  6,  the  scenario  is  essentially  the  same  except  for  additional 
branches  of  unstable  uniaxial  branches  that  exist  for  a  weak  external  field. 

Figure  3.2  depicts  the  stable  steady  states  along  with  some  nonparallel  unstable 
steady  states  at  two  selected  values  of  a  =  2.5  and  a  =  3.5.  The  metastable  steady 
states  become  available  at  about  iV>2.4  at  Q;  =  2.5,xo  =  l)M  =  0.1,i?  =  l  while  they 
exists  for  all  fV>0  at  a  =  3.5,xo  =  =  0.1,if  =  1.  The  stronger  dipole-dipole  poten¬ 

tial  facilitates  the  metastable  state.  We  note  that  the  value  of  the  order  parameter  in 
the  globally  stable  steady  state  remains  higher  than  that  of  the  metastable  state,  in¬ 
dicating  the  slight  loss  of  orientation  order  when  the  polarity  vector  is  aligned  against 
the  direction  of  the  external  field. 

4.  Conclusions 

We  have  extended  the  order  reduction  procedure  developed  in  [20]  to  complete 
the  study  of  steady  states  and  their  stability  in  extended  nematics  or  nematic  sus¬ 
pensions  when  the  dipole-dipole  and  the  excluded  volume  interaction  is  coupled  with 
the  external  fields.  The  main  results  are  summarized  in  the  following. 

•  In  the  absence  of  permanent  dipole  or  magnetic  moment,  both  the  polarity 
vector  and  the  external  field  are  shown  to  be  aligned  with  one  of  the  eigen¬ 
vectors  of  the  nematic  order  tensor.  As  a  result  of  this,  all  equilibria  can  be 
characterized  in  terms  of  three  scalar  order  parameters  and  a  rotational  group 
SO{2)  that  leaves  the  external  field  invariant.  Response  phase  diagrams  of 
the  order  parameters  with  respect  to  the  material  parameters  are  obtained 
together  with  the  stability  of  the  steady  states  determined  from  the  minimum 
of  the  free  energy  density.  Globally  stable  steady  states  are  uniaxial,  prolate, 
purely  nematic  or  polar  nematics  depending  on  the  strength  of  the  external 
field,  the  anisotropy,  the  strength  of  the  dipole-dipole  as  well  as  the  excluded 
volume  interaction.  A  family  of  biaxial,  metastable,  steady  states  may  exist 
with  the  major  director  aligned  in  the  direction  of  the  external  field  in  the 
range  of  higher  concentration.  The  polar  nematics  always  exist  in  pairs  with 
polar  order  ±si. 

•  In  the  presence  of  the  permanent  dipole  (or  magnetic  moment),  the  steady 
states  are  polar  nematics.  When  the  material  parameters  satisfy  the  condi¬ 
tions:  1^1  >XoE  and  3Nfj,<  \axoE\,  the  polarity  vector  and  the  external  field 
are  parallel  to  each  other  and  are  eigenvectors  of  the  nematic  order  tensor; 
while  the  conditions  are  not  satisfied,  nonparallel  steady  solutions  that  are 
in  plane  with  the  polarity  vector  and  the  external  field  can  exist,  but  are  un¬ 
stable.  In  any  case,  there  may  exist  up  to  two  stable  steady  states  which  are 
prolate,  uniaxial,  and  polar.  The  globally  stable  steady  state  has  its  polarity 
vector  parallel  to  the  external  field  direction  while  the  metastable  state  has 
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n=0.1  a=2.5  «=3.5 


Fig.  3.3.  The  polarity  vector  components  as  functions  of  N  for  a  =  2.5  and  0  =  3.5,  respectively. 
The  other  parameter  values  are  xo  =  solid  curves  are  stable  states.  The  stable 

and  metastable  steady  states  are  uniaxial,  prolate,  polar  nematics  with  the  polarity  vector  parallel  to 
the  external  field.  A  metastable  state  can  exist  at  higher  concentration  at  a  =  2.5  and  for  all  N  at 
0  =  3.5.  The  nonparallel  steady  states  in  which  the  polarity  is  not  parallel  to  the  external  field  can 
exist,  but  are  unstable.  The  nematic  order  of  the  globally  stable  steady  state  is  higher  than  that  of 
the  metastable  one. 
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its  polarity  vector  opposite  to  the  field.  The  former  exists  in  any  given  pa¬ 
rameter  regimes  while  the  latter  only  exists  at  regimes  of  high  strength  of  the 
dipole-dipole  interaction,  high  concentration,  and  strong  external  field.  The 
presence  of  the  permanent  dipole  breaks  the  symmetry  in  the  total  potential, 
yielding  the  two  stable  uniaxial  polar  nematics. 

This  study  completes  a  series  of  studies  on  the  potential  for  extended  nematics 
and  suspensions.  The  order  reduction  strategy  used  in  these  studies  works  for  any 
total  potential  given  as  a  finite  sum  of  lower  moments.  It  does  not,  however,  apply 
to  the  Onsager  potential,  where  the  potential  depends  on  infinitely  many  moments  or 
spherical  harmonics.  New  methods  would  have  to  be  invented  then. 

Appendix  A.  Free  energy  density  and  its  second  variation. 

A.l.  Free  energy  density  at  steady  states  in  order  parameters.  The 

free  energy  density  at  steady  states  in  order  parameters  when  the  permanent  dipole 
is  absent  is  given  by: 


A[f]  =  -kBT\nZ+^{asj+Nis^-sP  +  p'^)), 
Z=  f 


pasi  +  ^[(s— f  )(cos^6»— ^)  +  §sin^0cos2(/>] 


m. 


l|m||  =  l 


(A.l) 


Its  first  and  second  derivatives  are  given  by 

dA 


dsi 


=  aksT^si  —  {cos9)), 


r)A 

dA 


1-/3  3  2  A 

— -2(cos  d)j, 

^-^-^(sin^6»cosV)  )  , 


d^A 


=  akBT  —  a^ksT  {{cos^  9)  —  (cosd)^) 


dsl 

=  akBT  —  a^kBT 
d^A  SNakBT 


1  +  2S-I3 


■  -  Si 


dsids  2 

SNakBT 


{{cos^9)  —  (cos6i)(cos^0)) 

/  3a\  l-|-2s  — /3 

,cos  9) -Si - - - 


d'^A 

dsidp 


SNakBT 


( (cos  9  sin^  9  cos^  (j))  —  (cos  d)  (sin^  9  cos^  (j)) ) 


SNakBT 


(cosdsin^dcos^  (j))  —  si 


1-S  +  2P 


=  kBTN-  ((cos^  9)  -  (cos^  0)2) 


=  kBTN- 


4 

m'^kBT 


(cos^  9)  — 


{l  +  2s-(3f 
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d^A  ksTN 


dsdp  2  4 

keTN  m'^keT 


fp.  A 

4 


2  4 

9N^kBT,,_.  4,_4 


( (cos^  9  sin^  9  cos^  (\>)  —  (cos^  9)  (sin^  9  cos^  </>) ) 
cos^  0sin^0cos^(/))  — 


J.\  (l+2s-/3)(l-S  +  2/3) 


=  kBTN- 


( (sin"‘  9  cos"^  (/>)  —  (sin^  9  cos^  (j))  ^ ) 
9N^kBT  ''  o  „  ,  0^2 


(l-s  +  2/3)" 


[siw^cos"^  (/>)  —  ■ 


(A.2) 


At  the  isotropic  branch  (si  =  0,  s  =  0,/3  =  0),  the  second  variation  of  the  free  energy 
density  is 

(5^A|^^=o,s=o,/3=o  =  kBT[{a  -  a‘^/3)6^si  +  A^(l  -  N/5)d‘^s  -  N{1  -  N/5)ds6f3  + 
N{l-N/5)6^(3]. 

So,  the  isotropic  equilibrium  is  unstable  whenever  a  >  3  or  N  >5. 

For  the  nonzero  purely  uniaxial  nematic  equilibria,  s^0,/3  =  0, 

a(l  +  2s)\ 


6^A  =  kBTa  1- 


Ssf  +  kBTN{ds,  5/3)-C-iSs,  5(3)  ^ 


(A.3) 


where 

C  = 


(A.4) 


^(s)=  f  (A.5) 

Jo 

A.2.  Hessian  of  the  free  energy  for  dilute  polar  nematics.  In  the  dilute 
limit,  the  excluded  volume  potential  is  neglected.  The  Hessian  in  the  second  variation 
of  the  free  energy  density  in  steady  states  in  the  order  parameter  space  is  given  by 


H  =  a[(l-a(4±^-4?  -si))nn+(l  -  |(1  -s  +  2/3))nj_nx  + 
(1-  f  (l-s-/3))n*n*] 


(A.6) 


A.3.  Hessian  of  the  free  energy  density  for  the  concentrated  polar  ne¬ 
matics.  In  this  case,  both  the  dipole-dipole  and  excluded  volume  interaction  are 
included  in  the  potential.  We  denote 


M  =  (mm). 


A 


A 


12  = - ^[(ni)((Mii  —  M33),(M22  —  M33),Mi2,Mi3,M23)  — 

((ni(mf  —m§),ni(m2  —rn§),inrnim2,nirnim3,nim2m3)}j, 

(3JV)^  r/l  2  2  2  2  \T 

22  =  — ^ — [((mi  -m3, m2 -TO3, 77117712, mim3,m2m3) 

(mf  —1113,1712  — 1713,17111712,17111713,17121713))  — 

{Mil  —  M22,M22  —  M33,Mi2,Mi3,M23)'^ 

(Mil  —  M22,M22  —  M33,Mi2,Mi3,M23)]. 


(A.7) 
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Then,  the  Hessian  in  the  second  variation  is  given  by 

(a^((mm)-(m)(m))  Hi2\ 

(A.8) 

AJ2  A22 ) 

A. 4.  Proof  of  the  equivalence  of  the  hydrodynamic  equilibrium  and  the 
thermodynamic  equilibrium  in  this  model. 

Theorem  A.l.  The  steady  state  solution  of  the  Smoluchowski  equation  is  given  by 
the  Euler-Lagrange  equation  of  the  free  energy  density. 

Proof.  The  free  energy  density  for  homogeneous  flows  of  rigid  nematic  polymers 
or  rod/platelet  dispersions  is  given  by 


A[f]=kBT  [ 

d\\m\\  =  l 


W+'^f+vj 


dm, 


(A.9) 


where  Vi  and  44  are  the  intermolecular  potential  and  the  external  potential,  respec¬ 
tively.  The  Euler-Lagrange  equation  obtained  from  minimizing  the  free  energy  density 
functional  is 


/X  — A  =  0, 


(A.IO) 


where 

SA 

fj,=  —  =kBT[\nf  +  Vi  +  Ve]  (A. 11) 

is  the  extended  chemical  potential  and  A  is  the  Lagrange  multiplier  with  respect  to 
the  constraint 

I  /dm=l.  (A. 12) 

d|lm|l  =  l 

/X  =  A  is  the  governing  equation  for  the  thermodynamic  equilibria.  Clearly  it  yields 
solutions  of  the  Smoluchowski  equation. 

On  the  other  hand,  the  steady  state  solution  of  the  Smoluchowski  equation  is 
governed  by 


7^■(/7^/x)  =  0.  (A.13) 

We  multiply  the  equation  by  /x  and  then  integrate  over  the  unit  sphere  ||m||  =  l  to 
yield 

I  fiTZ- f'R.g,dni=  —  j  f\\TZfi\\'^dni  =  0.  (A. 14) 


This  implies 


/(m,^)7^/x  =  0. 


(A.15) 


Finally, 


TZfi  =  0. 


(A.16) 
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This  is  equivalent  to  /x  =  A.  □ 

Appendix  B.  The  stability  and  the  second  variation  of  the  free  energy 
density.  The  stability  of  the  steady  states  of  the  Smoluchowski  equation  in  the 
Doi-Hess  kinetic  theory  for  extended  nematics  is  studied  by  the  second  variation  of 
the  generalized  free  energy  density.  The  free  energy  density  for  the  material  system 
subject  to  an  external  field  is  given  by 


A[f]  =  kBT 


flnf+-V^[f]f  +  Vef  )dm 


(B.l) 


We  denote  by  fe  the  equilibrium  pdf  solution  and  Sfi  the  perturbation  of  the  equi¬ 
librium  solution.  Then,  the  second  variation  of  the  free  energy  density  is  given  by 

S^A[M=^-^kBT  [  fK  +  v,[f,]f,)dm.  (B.2) 

^  J\\in\\  =  l\Je  J 

Let  fi  =  fef 2-  The  second  variation  reduces  to 

(5"A=^fcBT(/|+f7[/e/2]/2)e,  (B.3) 

where  ((•))e  denotes  the  ensemble  average  taken  with  respect  to  the  equilibrium  pdf. 
Recall  that  the  equilibrium  solution  is  in  the  Maxwell-Boltzmann  form, 

=  (B,4) 

For  the  particular  intermolecular  potential  studied  in  this  paper,  the  intermolecular 
potential  can  be  rewritten  as 

Vi  =  -asim^  -  -  1/3)  -|-/3(m2  - 1/3)),  (B.5) 

where  M  =  s(nn  — 1/3) -|-/3(nxn^  — 1/3) -I-I/3  is  a  second  moment  tensor  of  the  pdf 
and  i  is  an  index  between  1  and  3.  Without  loss  of  generality,  we  set  i  =  l  in  this 
paper.  We  consider  the  perturbation 


f2=X-{x). 


(B.6) 


/2dm  =  0. 


(B.7) 


First,  we  focus  on  a  special  perturbation  that  perturbs  only  the  order  parameters: 

a;  =  simi  -I-  [siml  — 1/3)  -|-/3(m2  —  1/3)].  (B.8) 


The  equilibrium  pdf  is 

Z 


(B.9) 


Derive  the  Hessian  for  the  quadratic  term  with  respect  to  s,/3.  We  note  that 

s=  {ml -ml),  13=  {ml -ml). 


(B.IO) 
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We  parameterize  the  molecular  direction  m  in  the  frame  of  the  eigenvectors  of  the 
second  moment  tensor 


m=  (cos0,sindcos^,sindsin(/)). 

In  terms  of  the  angles, 

s  =  (cos^  6  —  sin^  9  sin^  4>))P=  (sin^  ^  cos  2(/()) . 
The  Hessian  of  the  free  energy  density  is  given  by 


d'^A 


3N 


=  2[(wf )  -  {m^y  -  ((mf )  +  —  (TOimm) :  (m^mm))], 

d'^A 


dsis 


d^A 

dsifi 


=  2{mi{ml  - 1/3))  -2{mi){ml  -  1/3)  -  {a{mi{ml  - 1/3)) (mf)  + 
3N 

—  l/3)mm) :  (mi mm)  +  {m'l){mi{m\  — 1/3))  + 

3N 

—^{miinm.) :  {(ml  —  l/3)mm), 

=  2{mi{m\  — 1/3))  —  2{mi){m\  —  1/3)  —  {a{mi{m\  —  1/3))  (ml)  + 
3N 

((m^  —  l/3)mm) :  (mi mm)  +  {ml){mi{m\  — 1/3))  + 


3N 


mimm) :  {{m\  —  l/3)mm), 


02  A 

—  =  2[{{ml  -  1/3)^)  -  {(ml  - 1/3))^  -  (a(m.(m?  -  1/3))^  + 
3iV 


"^^{{ml  —  l/3)mm) :  {{ml  —  l/3)mm))], 


d^A 


^  =  2((m?  -  l/3)(m^  - 1/3))  -  2((m?  -  1/3)) ((m^  - 1/3))  - 
asp 


(B.ll) 

(B.12) 


(B.13) 


{2a{mi{ml  —  l/3)){mi{ml  —  l/3))+3N{{ml  —  l/3)mm) :  {{ml  —  l/3)mm)), 


92  B 


^  =  2[((mi  -  1/3)^)  -  {{ml  - 1/3))^  -  {a{myml  -  1/3))^  + 
op^ 


3N 


{{ml  —  l/3)mm) :  {{ml  —  l/3)mm))]. 


Here  we  normalize  the  free  energy  density  by  ksT  and  i  =  l. 
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